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PERFORMANCE THEORY OF DIAGONAL CONDUCTING WALL MHD ACCELERATORS 

Ron J. Litchford* 

NASA - Marshall Space Flight Center 
Huntsville, AL USA 

The theoretical performance of diagonal conducting wall crossed field accelerators is m 

the basis of an infinite segmentation assumption using a cross-plane averaged generalized Ohm’s law 
for a partially ionized gas, including ion slip. The desired accelerator performance relationships are 
derived from the cross-plane averaged Ohm’s law by imposing appropriate configuration and loading 
constraints. A current dependent effective voltage drop model is also incorporated to account for 
cold-wall boundary layer effects including gasdynamic variations, discharge constriction, and 
electrode falls. Definition of dimensionless electric fields and current densities lead to the 
construction of graphical performance diagrams, which further illuminate the rudimentary behavior 
of crossed field accelerator operation. 


Nomenclature 

A cross-sectional area (= hw) 

A f slanted area (= A tan 0 ) 

b ion slip factor 

B magnetic field strength (induction) 

E electric field 

h channel height 

1 total current 

j current density 

k load factor 

L channel length 

M Mach number 

P power density 

R load resistance 

u velocity 

V voltage 

tv channel width 

p Hall parameter 

8 boundary layer thickness 

Y specific heat ratio 

A dimensionless effective voltage drop (Ej A uB ) 
tj electrical efficiency 

<p electric field direction (= tan 0 ) 
o scalar electrical conductivity 

2 effective electrical conductivity 

6 electric field angle 

0w diagonalization or wall angle (= n H - 0) 

£ dimensionless height parameter (= y/5) 

to conductivity ratio 

Q effective Hall parameter 


* Associate Fellow, AIAA. Copyright © 2003 by Ron J. 
Litchford. Published by the American Institute of 
Aeronautics and Astronautics, Inc. with permission. 
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Introduction 

Historically, interest in magnetohydrodynamic 
(MHD) devices has centered on their use as electrical 
generators in commercial central power plants and 
mobile burst power systems. The' primary attraction 
for central power is associated with the attainment of 
higher peak cycle temperatures, which point to 
significant improvements in overall plant efficiency. 
The attraction for mobile pulse power, particularly 
airborne military needs, derives from their intrinsic 
high power density characteristics. It should be 
recognized that equally important though less noted 
attractions exist for accelerator configurations, as well. 

Two major identifiable uses for crossed field MHD 
accelerators are as propulsive devices and hypersonic 
aerodynamic test facilities. Litchford et. al. discuss the 
operational attributes that are of particular significance 
to these applications and provide an in-depth historical 
perspective of their technological development. 111 

From a fundamental point erf view, the essential 
argument favoring utilization of a Loren tz force 
acceleration mechanism is the ability to avoid inherent 
physical limitations encountered with pure thermal 
approaches (viz, material thermal limits and 
ionization/dissociation losses). Simply put, it is more 
effective to transfer electrical energy into directed 
kinetic energy instead of first degrading it into thermal 
energy. Crossed field acceleration is of special interest 
in these cases because of its unique capacity for 
processing large amounts of power under conditions of 
high mass throughput. 

Past emphasis on generator configurations has led to 
a substantial theoretical and experimental basis for 
understanding their performance and operational 
nuances. However, a similar level of undo-standing 
has yet to be acquired for accelerator configurations, 
despite a strong growing interest in their potential 
aerospace applications. 
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The standard theoretical approach for describing the 
interaction of an electrically conducting gas with 
applied electric and magnetic fields relies cm 
application of a Cowling-Schl ter type generalized 
Ohm’s law for a partially ionized, electrically neutral 
gas. Powers et. al., for instance, used a cross plane 
averaged form of the generalized Ohm’s law and 
developed a graphical-based methodology for 
describing the general performance characteristics of 
both MHD accelerators and generators . 121 

This classical model has proven to be extremely 
useful in illuminating the basic operational behavior of 
MHD devices; but its practical utility is hampered by a 
critical limitation. That is, the theory assumes a 
constant near-electrode voltage drop, whereas 
experience has shown that the magnitude of the near- 
electrode voltage drop exhibits a significant load 
current dependency. 

Various theoretical refinements of this nature have 
been forthcoming, but solely within the context of 
generator configurations. Wu, for example, introduced 
an effective voltage drop parameter, accounted for its 
load current dependency, and examined the resulting 
impact on diagonal conducting wall (DCW) generator 
performance . 131 The MHD research group at the Max- 
Plank-Institut f r Plasmaphysik made a similar attack 
on the problem using a cross plane averaged Ohm’s 
law, as specialized for a DCW generator 
configuration . 141 


It is widely recognized, of course, that the current 
and electric field structure in MHD devices is 
inherently three-dimensional and spatial separation of 
physical processes is not warranted in general. As a 
matter of recourse, one could mount a brute force 
attack cm the problem and perform a complete three- 
dimensional numerical analysis of the flow and 
electrical structure in the duct, assuming adequate 
computing resources are available. Alternatively, one 
could invoke an infinite segmentation assumption (i.e., 
the stream wise variation in electrical properties is 
small in comparison to transverse variations) and 
couple an approximate cross plane electrical model 
with a three-dimensional flow analysis, as exemplified 
for generator configurations by Bityurin et. al . 151 and 
Ahluwalia et. al . 161 Each approach has it merits and 
drawbacks, and experience has shown that each level 
of approximation is useful in fulfilling certain design 
and analysis needs. 

The purpose of this paper is to reexamine the 
classical theory and extend it for cross-field accelerator 
configurations with inclusion of a current dependent 
near-electrode voltage drop model. This particular 
refinement of the classical theory, while generally 
recognized as a straightforward theoretical extension, 
has not explicitly appeared in the literature, and it is 
believed that the resulting development can yield 
practical insight into the basic operational 
characteristics of these important devices. Moreover, it 
is hoped that this work may serve as a convenient and 
compact resource for future design practitioners. 
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(a) Linear Hall Accelerator 



(c) Series Connected Diagonal Accelerator 


(b) Segmented Faraday Accelerator 



(d) Diagonal Conducting Wall Accelerator 


Figure 1: Composite graphic of alternative cross-field MHD accelerator configurations. 
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Crossed Field Accelerator Configurations 

Various linear crossed-field MHD accelerator 
configurations are conceivable, as illustrated in Fig. 1, 
although the optimal configuration depends upon the 
ultimate application need. The Hall configuration, for 
example, is generally more effective for low-density 
flows (i.e., very high Hall parameter condition) 
whereas the segmented Faraday configuration yields 
superior performance for high-density flows. A 
significant disadvantage of the Faraday configuration, 
however, is the separate power conditioning required 
for each electrode pair leading to a heavy, complex, 
and expensive system. 

Alternative two-terminal loading schemes have been 
proposed to avoid the complications of multi-terminal 
connections while also attempting to compensate for 
the effects associated with a finite Hall parameter. De 
Montardy, for instance, suggested the series connected 
scheme in which a segmented Faraday channel is 
externally diagonalized. 171 Dicks later extended this 
approach to a diagonal conducting wall (DCW) 
configuration in which slanted window-frame-like 
electrode elements are stacked with thin insulators to 
form a complete channel. 1 ®' 91 The DCW configuration 
not only simplifies fabrication and improves strength 
but provides superior performance to the series 
connected device by allowing current to flow to the 
sidewalls. 131 

Hence, theoretical development is focused on the 
DCW accelerator configuration, since it has the 
greatest practical relevance. The performance 
characteristics of the alternative crossed field 
configurations can be subsequently deduced as special 
or extreme cases of the DCW configuration. 

Generalized Ohm’s Law 

The local current density j in a partially ionized, 
electrically neutral gas may be determined from a 
Cowling-Schl ter type generalized Ohm’s law. 110-141 
For application to typical MHD devices, it is useful to 
adopt the simplified form derived by Rosa which 
includes ion slip but neglects electron pressure 
gradients: 1131 

j = oE'-— jxB+-^^-(jxB)xB (1) 

B B 

Here, a is the scalar electrical conductivity, E ' is the 
generalized electric field intensity relative to axes 
moving with gas velocity u, B is the magnetic field 
strength, and 0, and $ are the electron and ion Hall 
parameters, respectively. The Hall effect and ion slip 
are introduced by the second and third terms on the 
right hand side. 


As a practical matter, we define the generalized 
electric field E' in toms of the electric field in the 
laboratory frame E, the motional emf uxB, and an 
equivalent electric field Ej which accounts for the 
near-electrode voltage drop V d (i.e., E d = V d /h) inside 
the boundary layer: 

E'=E+ uxB-Ej (2) 

Note that the equivalent electric field is applied only in 
a direction transverse to the magnetic field vector and 
the stream wise velocity vector. It is always directed in 
opposition to the transverse current, j y . 

For the special case in which no current flows along 
the magnetic field lines (j z = 0), it is readily shown that 
Eq. (1) reduces to the compact form: 

j = IE'-— jxB (3) 

B 


where we have introduced two new parameters 
a a 
l+A/O+b 


z= 


(4) 


Q = 


B 0 

i +AA 


(5) 


In many cases, the ion slip factor b = 0 e 0 is negligible 
and we obtain further simplification as £-*tr and 


Mean Ohm’s Law 

Many of the interesting phenomena that occur in 
MHD devices such as velocity overshoots in the 
sidewall boundary layer, flow asymmetries, and 
generation of intense secondary flows in the cross 
plane are three dimensional in nature. Therefore, high 
fidelity performance predictions ultimately require 
three-dimensional numerical analyses of the combined 
flow and electrical structure. On the other hand, 
simplified approaches based on spatial decoupling or 
reduced spatial dimensions can also yield useful results 
if all critical physical effects are accounted for in an 
appropriate fashion. 

Exact theoretical treatment requires consideration of 
the fundamental equations of MHD and Maxwell’s 
electromagnetic equations, thereby leading to a system 
of coupled differential equations. Here, we invoke the 
infinite segmentation assumption and assume that all 
parameters in the generalized Ohm’s law can be 
effectively represented by their cross plane averaged 
values and are dependent on the x-coordinate, only. 

The dimensionless voltage drop parameter A = E d 
/uB = V d /uBh, which incorporates all voltage losses 
associated with the electrode boundary layer, may be 
introduced into the generalized Ohm’s law to obtain 
the following component relations: 
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l=XE x -nj y ( 6 ) 

j y ='LE y -I l uB(l+A)+Slj x (7) 

Equations (6) and (7) may also be combined to obtain 
the inverse expressions for j x and j y in terms of the 
local electric field components: 

(8) 

Performance Model 

Diagonal Loading Constraint 

Completion of the electrical model requires 
consideration of external loading conditions. In the 
case of a diagonally linked device, as illustrated in Fig. 
2, the electric field is forced to align perpendicular to 
the diagonal short such that 

E y /E x -Xan9 = q> (10) 

where 0 is the electric field angle and q> is the electric 
field direction. It follows that the diagonalization or 
wall angle is given by 0, = it 12 - 0. The sign of <p 
depends on the mode of operation. For an accelerator, 
uxB opposes jy, 0 < it 12 < 0, and <p> 0. For a 
generator, uxB aligns with jy, it 12 <0<it, and q>< 0. 

The two terminal load current / for an MHD device 
with diagonally linked electrode pairs is given by 

7= f jdA,=j nA, (11) 

JA f 

where the integration is over the entire slanted area A f . 



Figure 2: Orientation of field vectors in a 

diagonally connected cross-field accelerator. 

In component form: 

I = (j x + tan0j y )A = (j x + <pj y )A (12) 

Electrical Parameter Relationships 

A complete set of equations now exists for 
determining the cross plane electrical characteristics. 
In general, either / or E x may be specified allowing the 
remaining unknown electrical parameters to be 
deduced from Eqs. (6), (7), (10), and (12) in 
conjunction with appropriate material functions. 

The resulting performance relationships fix' DCW 
accelerators are summarized in Table 1 using both E x 
and 7 as the independent variable. The governing 
relations for generator performance are recovered when 
both q> and A are less than zero and 7 = -EJJR, where 
R is the load resistance and L is the channel length 
over which the load is applied. Note also that the 
governing relations for a linear Hall channel 


Table 1: Summary of Electrical Parameter Relationships for Diagonalized Cross-Field MHD Devices. * 


Electrical Parameter 

In terms of applied electric field, E x 

In terms of applied electric current, I 

Jx 

(l-n<p)lE x +OLuB(l + A) 

(l-n<p)l/A+<pE.uB(l + A) 

i+n 2 

X + q? 

jy 

(n+V»)ZE,-ZuB(l+A) 

1+n 2 

(n+9>)7/A-IuB(l+A) 

E x 

E x 

(l+£2 2 )7/A+(p-n)l«fl(l + A) 
E(lV) 

Ey 

<pE x 

q>E x 

I 

(l+<p 2 )l£ x A+(n-<p)?.uB(l+A)A 

l 


1+fl 2 


* The sign of <p and A depends on the mode of operation. Generator Mode: <p< 0 and -7 <A <0. 
Accelerator Mode: (p > 0 and A > 0. 
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configuration are recovered in the extreme case <p=0. 
Mean Power Parameters 

Ultimately, we seek to deduce suitably averaged 
power parameters from the resulting cross plane 
electrical model. The quantities j-E and (jxB) x , for 
instance, may be used to determine the power density, 
push power, dissipated power, and electrical efficiency 
of a device. 

Power Density 

The electrical power density at any cross section of 
the duct is defined by 

P = l-E = -!-jEl + E) (13) 

A f 

Eliminating E y using Eq. (10) yields 

P = -LJe] + <f? E] = (14) 

A f A f 

But Af = A^Jl + <p 2 and we obtain the working form: 

P = ( 15 ) 

A 


Push Power 

The stream wise Lorentz body force component at 
any cross section is defined as 

(J*B) x = j y B (16) 

And the push power associated with this Lorentz body 
force is given by 

P p = u-(jxB) x = uj y B (17) 


Dissipated Power 

The power density of an MHD device must exceed 
the Lorentz force push power since the internal 
resistance to current flow leads to ohmic heating. This 
can be demonstrated by taking the dot product of j with 
the generalized Ohm’s law, as defined by Eq. (3), 


j-j = Z [j-E+j.(uxB)-j-E,]-£j.(jxB) (18) 

B 

Because j z = 0, the last term is identically zero and we 
obtain 


= j E-u jxB- j E, = P-P p -j y uBA 


09 ) 


where j-j/Z is the power dissipated in the core flow 
region and j-E rf = jyE d = j y uBA represents the power 
dissipated in the electrical boundary layers. Therefore, 
the power dissipated over the entire cross plane is 
given by 

/> =ki+j yU BA = P-P r (20) 


Electrical Efficiency 

The electrical efficiency of an accelerator is simply 
the ratio of the push power to the applied power 
P f u(jxB), uj y B 

* P j-E IEjA 


( 21 ) 


It follows that the generator efficiency is defined as the 
reciprocal of the accelerator efficiency (ij g = l/ij a ). 


Power Relationships 

By combining the results in Table 1 with the 
relationships above, it is possible to express the mean 
power parameters in terms of the independent variables 
E x or /. These results are summarized in Table 2. 


Generalized Vector Diagram 

It is generally instructive to examine the qualitative 
behavior of the internal electric fields and current 
densities undo* the assumption that MHD effects have 
only a slight effect on flow properties. In this case, the 
generalized Ohm’s law can be used to construct 
representative vector diagrams for diagonalized 
devices, as shown in Fig. 3 for fixed values of <p. 


Table 2: Summary of Mean Power Relationships for Diagonalized Cross-Field MHD Devices. * 


Power Parameter 


P = j-E 
Pp = u-(jxB)^ 


In terms of applied electric field, E x In terms of applied electric current, / 


(l+^)+(£2-^)uB(l+A)/E, 


(l+Q 2 )+(^-Q)Z«B(l + A)iV^ 

1 + Q 2 

IA 2 

1+p 2 


uBl\{n+<p)-?MB{U-A)All 
A l + p 2 


I£uB 


■(n+ 


1+Q 2 


„ «(jxB), 1 uB 

[Sl+<p)-uB(\+A)/E x 

ZuB 

(a+<p)-ZuB(l+A)A/l 


j-E n s E x 

(l ^ + <p 2 ) + (Q- <p)uB (1 + A )/E x 

I/A 

(l+Q 2 )+(p-Q)Z«£(l + A)A/7 



* The sign of 0 >and A depends on the mode of operation. Generator Mode: <p < 0 and -1 <A < 0 . 
Accelerator Mode: q> > 0 and A £0. 
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(d) open circuit accelerator 


(e) zero efficiency accelerator 


(f) ^neutralized accelerator 


Figure 3: Vector diagrams for crossed field diagonal MHD devices. Each diagram depicts the relative 
orientation of the current density vector and electric field vectors based on the generalized Ohm's law 

^+q(-^x— j = E' and the generalized electric field relation E' = E+uxB-E rf . 

£ B J 


To construct these diagrams, the generalized Ohm’s 
law defined by Eq. (3) is put into the form 



where we immediately recognize that the two terms on 
left hand side must be orthogonal. Thus, each of these 
terms represents one side of a right triangle with the 
hypotenuse E' inclined at an angle of tan' 1 Q. The 
generalized electric field must simultaneously satisfy 
the vector relation defined by Eq. (2). 

The top three diagrams illustrate the load current 
dependence of a diagonalized generator ( <p< 0). In the 
open circuit case (a), there is no current flow to the 
load, and the internal current density vector must align 
with the diagonal link, which defines planes of 
constant electric potential. Thus, the current is closed 
through the diagonal short, and Ey exactly compensates 
for the combined uxB and E</ induced potentials. Note 


that j z is directed upstream under these conditions. As 
the load current is increased, the current density vector 
rotates out of the plane of the diagonal link until it 
aligns with the uxB vector. At this point, j x = 0 and 
we arrive at the Hall current neutralized condition, (b). 
If the load current is further increased, the current 
vector continues its rotation and the direction of Hall 
current flow is reversed. Ultimately, the current 
density vector becomes perpendicular to the plane of 
the diagonal link as the electric field goes to zero at the 
short circuit condition (c). 

The bottom three diagrams illustrate the load 
current dependence of a diagonalized accelerator {<p> 
0 ). When the applied axial electric field is zero and 
there is no load current, we obtain the open circuit 
accelerator loading condition (d). The internal current 
density vector corresponds to the generator shot 
circuit case, under these unique circumstances. Note 
that j y < 0 (braking regime) and that the current flow is 
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closed through the diagonal short. With increasing 
applied axial electric field, the current density vector 
rotates out of the plane of the diagonal link. The 
device remains in a braking regime, however, until the 
applied field becomes large enough to force j y to zero, 
at which point we attain the zero efficiency accelerator 
condition (e). Further increases in applied current 
leads to continued rotation of the current density vector 
and increasing values of j y . Eventually, one reaches a 
j x neutralized accelerator condition (f). Increasing the 
applied current beyond this point leads to a reversal in 
the direction of axial current flow while the magnitude 
of the current density vector continues to grow. 


Current Dependent Voltage Drop 


The previously developed performance relationships 
hold for any fixed value of A. However, extensive 
experimental research with MHD generators has 
clearly shown that the near-electrode voltage drop is 
linearly dependent on the transverse current density. 116 * 
181 That is, any increase in load current leads to a 
proportionate increase in voltage drop. This behavior 
is assumed to extend to accelerator operation, as well. 

The linear relationship between Vj and j y permits 
the construction of a simple two-parameter model 

W 




I'- I' 


(i-o+v; 


(23) 


where the prime and double-prime superscripts 
indicate two distinct loading conditions. It is generally 
more convenient to express this relation directly in 
terms of A 

Ta'-a^ 


A = 


I' -I' 


(/-0+ A' 


(24) 


Selection of appropriate reference load conditions 
depends on the mode of operation. For a generator, it 
is most natural to adopt the open circuit (A' = A „) and 
short circuit (A" = A,) conditions on the extreme ends 
of the operational load line. For an accelerator, it is 
convenient to adopt the zero power condition, which 
corresponds to a short circuit generator (A' = A z = A,), 
and the j y neutralized condition (A" = A, = 0). The 
selected reference conditions are summarized in Table 
3. The implications of this model are now examined 
for both generator and accelerator modes. 


Table 3: Reference conditions for variable A. 


Mode 

A' A' /' 

r 

Generator 

> 

> 

II 

o 

h 

Accelerator 

A Z = A , A l> = 0 7 Z 

7 „ 


Generator Mode 

At the short circuit generator condition, E x = 0 , f' = 
/„ and A" = A,. From Table 1, we therefore deduce 
r = I l + A,)(fl-g) 


l +£2 


(25) 


Elimination of 7" using Eq. (24) immediately yields an 
expression for A in terms of / for a generator 

'(l + tf)(A.-A.) 


A = A - 


(Q-p)(l+A,) 


ImBA 


(26) 


This relationship may be used to eliminate A in the 
performance relations of Tables 1 and 2. 


Accelerator Mode 

The zero power open circuit accelerator condition 
corresponds to the short circuit generator condition 
with 9 > 0. In this case, E x = 0, T - l z , A' = A z = A, 
and we deduce that 


( 27) 

l^fl 

The point of transverse current density neutralization 
corresponds to a zero efficient accelerator where j y = 0, 
I" = /„, and A" = A„ = 0. This implies 

ImBA 


/'=/=- 


tl+<p 


(28) 


Eliminating T and I" in Eq. (24) gives an expression 
for A in terms of / for an accelerator 


A = A. 


1- 


(Q + rt(i + a’)^-(n ! -y’)a + 4,) 

(i+n 2 )-(n i -9 2 )0+A > ) 


(29) 


This relationship may be used to eliminate A in the 
performance relations of Tables 1 and 2. 


Performance Diagrams 


It is possible to construct an informative 
performance map for DCW MHD devices following 
the general methodology first outlined by Powers et. 
al., who assumed a constant voltage drop in the 
electrode sheath layer. 121 The significant modifications 
associated with a variable effective voltage drop model 
were later noted and described by Wu for the DCW 
generator case. (3] Here, Wu’s analysis is explicitly 
extended for DCW accelerator operation, as well. 

Dimensionless Performance Parameters 


As a first step, the governing performance 
relationships are simplified by defining the following 
dimensionless quantities: 

j „ E „ / _ P 


y- J c _ JL t — * V = ■ 

J luB ’ uB ’ ImBA ’ Iu 2 B 


2 r>2 


( 30 ) 
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Substitution into the component form of the 
generalized Ohm’s law then yields 

J,=I^r[£,+^,+n (l +A)] (31) 

+ (32) 

It also follows that the electric field components can be 
, stated explicitly in terms of current density 

£ X =J X +QJ, (33) 

£ f =-aj x +J,+( 1+A) (34) 

An expression for the dimensionless load current is 
obtained from substitution of Eq. (30) into Eq. (12) 

1 = J,+<PJ, (35) 

And the dimensionless laboratory power density 

follows directly from substitution of Eq. (30) into Eq. 
(13) 

V = J£ = J x £ x + J,£ y (36) 


Eliminating £x and £y using Eqs. (33) and (34) yields 
the more convenient form 


V = J 2 x +J 2 +(l+A)J y (37) 


Furthermore, we may define the dimensionless push 
power as 


V = - p - 
” I u 2 B 2 


_ UjyB jy 


'Lu 2 B 2 'LuB (38) 

And the electrical efficiency may be expressed in toms 
of dimensionless parameters by combining Eqs. (37) 
and (38) 



£ 




(39) 


It is readily apparent that the relationships in Tables 1 
and 2 can all be put in dimensionless form using Eq. 
(30); the results are obvious and are not tabulated. The 
current dependent effective voltage drop effect is 
accounted for by using either Eq. (26) or (29) in 
dimensionless form. 

As a prelude to the construction of the graphical 
performance map, it is instructive to examine 
performance parameter variations as a function of the 
dimensionless axial electric field, £ x . Figure 4, for 


example, summarizes calculated parameter variations 
assuming representative values for electrical/flow 
conditions, material functions, and wall angle. 
Starting from the open circuit generator condition and 
gradually reducing the load impedance, we observe the 
sequential occurrence of all loading conditions 
depicted in Fig. 3. 

For an open circuit generator, T = 0 and Eq. (35) 


yields the constraint = - <pj y . 0 * We therefore infer 


8 



Figure 4: Variation of performance parameters 

based on representative values for flow/field 
conditions , material functions ; and wall angle . 
Assumes <p=±l t B=4T y £ - 100, S2 = 2, u ~ 3 
km/s y Ao = -OJy and A* = -0.02. 

that ffx.o < 0 since 9 < 0 and J ytC < 0 for the assumed 

magnetic field direction. In this case, the internal 
current is forced parallel to the wall angle and is 
completely shorted through the diagonal linkage. As 
the load impedance is decreased below a critical level, 
however, the total current becomes nonzero and the 
negatively directed axial electric field begins to fall in 
magnitude. 

If the load impedance is decreased all the way to 
zero, we arrive at the short circuit generator condition 
corresponding to £, = 0. Here, application of Eq. (33) 

yields the constraint = - f 1J y ^ which implies J x>5 
> 0 since Q > 0 and J y , s < 0. In this case, the internal 

current flow is forced perpendicular to the wall angle, 
and no current flows through the diagonal linkage. 
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Consideration of these open and short circuit 
constraints clearly implies the existence of a Hall 
current neutralized condition (i.e., J x = 0) on the 

generator load line. It is also clear that maximum 
power extraction and electrical efficiency must also 
occur at some point between open and short circuit 
conditions. This point may or may not coincide with 
the Hall current neutralized condition. Note that the 
magnitude of A is largest at open circuit loading where 
the magnitude of J y is greatest. 

Simply reversing the wall angle (9 > 0) at the 
generator short cirpuit condition yields the open circuit 
accelerator condition. If power is then gradually 
applied such that a positive axial electric field is 
impressed on the accelerator, we enter a regime where 
the transverse voltage is insufficient to overcome the 
induced uxB potential, and J y remains negatively 
directed. This is commonly referred to as the braking 
regime where all applied power goes into joule heating 
of the working fluid. Eventually, the applied voltage 
becomes high enough to neutralize the transverse 
current at which point A -> 0. 

Increasing the applied electric field beyond this 
point yields positive push work and flow acceleration. 
The transverse current and total current grow in 
magnitude, more power is delivered to the device, and 
A steadily increases. At the same time, axial current 
flow gradually decreases until the J x neutralization 
condition is reached. Additional increases in applied 
field and power beyond this point simply reverses the 
direction of J x . It is important to note that optimal 
accelerator electrical efficiency occurs at an applied 
electric field much less than that required to neutralize 
axial current fiow. 

Performance Map 

The operational attributes of MHD devices can be 
exhibited in the most vivid manner on a performance 
map in the plane of dimensionless current density and 
electric field. Following the methodology of Powers et. 
al., m as modified by Wu for current dependent voltage 
drop, 131 we complete the square in Eq. (37) to obtain 
the family of circles defined by 

Jl+[j y +(l+A)/2j =7>+[(l+A)/2] 2 (40) 

These circles are centered at J x = 0, J y = (l+A)/2 with 
a radius of [P+ [(1+A)/2] 2 } 1/2 . 

Transformation to the £„ £y axes is defined by Eqs. 
(33) and (34), and the origin of the 6„ £, axes in the 
J„ J y plane follows directly from Eqs. (31) and (32) 


Q(l+A) 
”'° 1+Q 2 


j = _(L^1 (41) 

l+Q 2 () 


The locus of points defining this origin can be 
expressed independent of Q by combining Eqs. (33) 
and (34) under the condition £* = £>, = 0 

3 14-0 + [^-o + (l+A)/2] 2 =[(l+A)/2] 2 (42) 


Note that the circles defined by Eq. (40) reduce to 
the loci defined by Eq. (42) at zero power density. 
Thus, the generator boundary coincides with the locus 
of the origin of the £„ £y axes, and all power 


generating circles lie inside this boundary since V < 0 
for a generator. Further insight may be gained through 
examination of extreme generator loading conditions. 

At generator short circuit, the Hall field goes to zero 
and Eq. (33) gives the constraint 

J x ,=-OJ y „ => tm l (-]/£!) = jjj„ (43) 

Because Q > 0 and J y < 0, given the assumed magnetic 
field orientation, the short circuit operating domain 
resides in the fourth quadrant of the current density 
plane. 

At generator open circuit conditions, the total 
current falls to zero and Eq. (35) gives the constraint 
J,, =-<PJ y ,o => tan- 1 (-1 l<p) = Jy,lJ x , (44) 

Because 9 < 0 and J y < 0 for the generator mode, the 
open circuit operating domain resides in the third 
quadrant of the current density plane. 

Clearly, the power density vanishes at both open and 
short circuit conditions, but the critical discriminating 
factor to note is that A varies with current over the 
entire load range (-1 < A* < A, < 0). Therefore, the 
radius of the zero power circle depends on the load. 
The zero power loci for short and open circuit loading 

follows from Eq. (40) with V=0 

Jl +[j y , +(l+A,)/2j =[(l+A,)/2] 2 (45) 

+(l+A,)/2] 2 =[(1+ Aj/2] 2 (46) 

Using these expressions, it is possible to construct 
the short and open circuit half-circles. The short 
circuit half-circle is centered at J x = 0, J y = (l+A,)/2 
with radius (l+A,)/2, and the open circuit half-circle is 
centered at J x = 0, J y = (l+AJ/2 with radius (l+AJ/2. 
These loci are illustrated in Fig. 5 where the load line 
connects in accordance with the constraints defined by 
Eqs. (43) and (44). Because A„ < A„ the radius of the 
open circuit half-circle is smaller than the radius of the 
short circuit half-circle and the center points are 
shifted. 
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Figure 5: Performance map for diagonal conducting 
wall MHD devices with current dependent voltage 
drop effect. 

The origin of the £„ £y axes falls on the short circuit 
half-circle at the points defined by Eq. (41). Inspection 
of Eq. (34) indicates that the £, axes passes through the 
origin of the J„ J y axes since A = 0 at this condition. 
Various points of interest are indicated cm the diagram. 

Combining the DCW electric field constraint 
defined by Eq. (10) with Eqs. (30), (33) and (34) yields 
an expression for the operating load line in the form 

(Cl+q>)J x +(to<p-\)J y =1+A (47) 

It is now possible to eliminate A using either Eq. (26) 
for a generator or Eq. (29) for an accelerator. In both 
case, one ends up with an equation for a line in the J„ 
J y plane since A is linearly related to J x and J y . The 
resulting operating lines are illustrated in Fig. S. 

In the generator case, the operating line passes 
through the short and open circuit load conditions, as 
determined by Eqs. (43) and (44), respectively. In the 
accelerator case, the operating line originates from the 
generator short circuit loading point on the zero power 
half-circle. This line also projects through a point cm 
the half-circle defined by the value of q>. 

The V circles for accelerator operation are obtained 

from Eq. (40) using Eq. (29) to eliminate A. Note that 
the accelerator operates in the dissipative braking 
regime between open circuit and zero efficiency 


10 


conditions. Thus, positive push work is obtained only 
when the power density exceeds that associated with 
the zero efficiency condition where J y = q a = A„ = 0 



Portions of the constant-? 7 circles lying within the 
positive push work regime are illustrated in Fig. 5. 
Inspection of Eq. (40) reveals that the radii of these 
circles are larger than would occur for an ideal device 
in which there is no parasitic voltage drop (A = 0). 
Thus, adolitiomal power must be applied in a real device 
to achieve a desired intensity of acceleration. 

It is important to note that the optimum accelerator 
efficiency occurs with a finite axial current. Increasing 
the accelerator power density until J x vanishes will 
inevitably reduce the electrical efficiency below the 
optimal value. In fact, to accelerate plasma efficiently 
(with little beat production), it is desirable to maintain 
gentle acceleration levels by keeping the back e.m.f. 
only slightly less than the applied electric field. In 
practice, the added length and weight associated with 
gentle acceleration must be traded against the electrical 
inefficiencies encountered with high push power. 

For a generator, the maximum power condition can 
differ considerably from the ma ximum efficiency 
condition. Furthermore, maximum generator 
efficiency operation also occurs with a finite axial 
current, as observed in Fig. 4. It should be pointed out 
that achieving optimal power generation for a given 
value of Hall parameter depends on careful selection of 
the wall angle. This is examined for the sake of 
completeness. 

By definition, the olimensionless power density may 
be expressed in the form 

If we use Eq. (26) to eliminate I in the appropriate 


expressions for J„ Jy, and 6„ we o)btain the 
dimensionless generator power density in terms of A 


(Q-p 2 )(l+A,) 2 j 

a-a„ ) 

’a- A. 1+a' 

( 1+ ^ 2 )( 1+£ ^ 2 ) 

Ia-a J 

A -A 1+A 

L * O s J 


(50) 


The maximum power density is defined by the 
constraint dV/dA = 0 whit* yields the relation 


Am* =(A 0 +A,)/2 (51) 

Substituting this result into Eq. (50) then yielols 
-(n-f*) (1+A,)(1+A,) 

“ (l+^Xl+O 2 ) 4 
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Thus, V max depends cm <p, and the optimal wall angle 
for a given Hall parameter is determined by the 
constraint dVmJd 9 = 0 which leads to the design 
criterion 

Cl<p = - 1 (53) 

If this criterion is enforced, we find the maximum 
possible generator power density to be 

^ (n ,— ,)=-(l+A,)(l+Aj/4 (54) 

Electrode Boundary Layer Effects 

The sheath layer voltage drop arises from cold wall 
boundary layer effects including gasdynamic 
variations, discharge constriction, and electrode falls. 
Because fluid temperature and velocity vary rapidly 
approaching the wall, the conductivity and induced 
voltage also exhibit strong variations in the thermal 
and momentum boundary layers, respectively. 
Furthermore, the rapid decrease in temperature near 
the electrode surface leads to breakdown and a sudden 
switch from diffuse current transport to constricted 
arcs. Current attachment to cold electrodes almost 
always occurs through a thin layer of short arcs. 

The basic boundary layer features are depicted in 
Fig. 6 for an accelerator configuration. For simplicity, 
the thermal and momentum boundary layers are both 
assumed to have a thickness 8. In addition, current 
constriction is indicated at the electrode surface where 
attachment occurs through short arcs spanning the 
anode thickness 5 a and the cathode thickness 5 C . Note 
that the illustration shows a concentration of current 
due to the Hall effect at the upstream end of the 
cathode and the downstream end of the anode in 
accordance with experimental observations for 
accelerators.' 19,201 

The effective electrode voltage drop V d is defined as 
the potential difference by which the voltage between 
opposite electrodes is increased (due to boundary layer 
effects) with respect to the folly homogeneous case. 
The increase in V d is quantified through an integration 
of the transverse electric field defect across the 
electrode boundary layers. For convenience, V d is split 
into the following separate components: 

V d =V m +V e +V, (55) 

where the anode voltage fall V a , cathode voltage fall V c , 
and gasdynamic voltage drop V, are defined by 

*£*(*>-**)* (56 > 



Figure 6: Illustration of current transport and 

boundary layer dimensions for a diagonal conducting 
wall cross-field accelerator. 


Note that the 0 subscript refers to the homogeneous 
core flow conditions. 

Factors contributing to the anode and cathode 
voltage falls include the arc attachment region and the 
potential associated with electron emission from cold 
electrode material. In general, the voltage drop in the 
arc columns is negligible since the arc length is small 
compared to the boundary layer thickness, and the 
anode electrode fall is negligible since it doesn’t have 
to emit electrons. Thus, V a * 0. The cathode does emit 
electrons, however, and the voltage drop associated 
with this process can only be approximated in an 
empirical fashion. For copper, V c » 20 volts. 

The gasdynamic voltage drop may be evaluated by 
applying the generalized Ohm’s law in the boundary 
layer to eliminate E y . In this particular case, the 
appropriate form for E y may be obtained from Eq. (9) 
by neglecting A 

E t =^j^j y -ClE I+ uB (57) 


Assuming Q is invariant, V g takes the form 

v *-£fe£ =rt-1 - r - A 


[z/l+Q; 
* [[^1+ 




w g 

■H 

H 


/ ^ + M o® 

Zo/l+Oo ° 


jy 

So/l+flo 


1‘ 


(58) 


And it is convenient to split V g into two parts 


II 

if-: 

jrfy 

(59) 

V u =2u 0 B( 


ijrfy 

(60) 

J 

H«o 
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where Vi and V u account separately for the decrease of 
X and induced voltage in the boundary layer and X' = 
Z/d+flo 2 ). These integrals may be further simplified 


by defining % = y/8 and assuming 8 a = 8 C 


v = 2Sjy f 

;[#-] 

k 

(61) 

V u =2Su 0 Bj 

■U-i 

f<\u n 

dt 

(62) 

The lower limit of integration ^ 

= 8 is the distance 


from the electrode where current becomes constricted. 

Evaluation of these integrals is facilitated by the 
introduction of model boundary layer profiles for 
velocity and temperature. 1211 For example, the velocity 
profile in a turbulent boundary layer is well 
approximated by the power law u/uo = (y/8) 11 ” where n 
* 7. If this relation is substituted into the integral for 
V u and - 0 is used as the lower limit of integration 
(since u goes to zero at the wall), we obtain 

«* 2u 0 B6/(l+n) (63) 


Evaluation of the integral for Vi is complicated by 
the fact that X'o/X' increases steeply near the cold 
electrode surface, which makes the integral extremely 
sensitive to the lower limit of integration. It is 
assumed that electrical breakdown leading to arc 
attachment occurs when j/L' exceeds the critical 
electric field strength E c . In this case, the limit ^ 
corresponds to a critical conductivity ratio cot 


G>c 


K J,/E, 


y' 


L' 


(64) 


Experience shows that the turbulent thermal 
boundary layer may be approximated by the power law 
T/T 0 = (y/S) llm where m = 18 - 25. Furthermore, it has 
also been established that the electrical conductivity 
can be approximated by the simple power law X7X' 0 = 
(T/T 0 f(p/po) b where a - 10 - 20 and b = -1/2. Thus, 
for negligible pressure variations, the integral for V £ 
can be expressed in the simplified form 

V'-Hhttr 1 -- a* <«> 


The calculation procedure is as follows. For a given 
value of jyfL'o and a specified E„ it is possible to 
compute at which may be used in turn to compute 
TJT 0 = cot 1/a . Then, the lower limit of integration may 
be determined from the relation ^ = (TJT 0 ) m and the 
integral of Eq. (65) may be evaluated. For illustrative 
purposes, calculated values for ^ and Vj/28 are shown 
in Fig. 7 as a function of jJL’o assuming E c = 40 
kV/cm, n = 7, m = 18, and a = 10. Note the linear 
dependence of Vj on j y , in accordance with empirical 
observations. 


12 



Figure 7: Variation of%c and Vi/28 as a function 
of j/L' o . Assumes E e = 40 kV/cm, n = 7, m - 18, 
and a = 10. 


Practical Design Constraints 

In a real device, practical design constraints limit 
the available operating range. The axial electric field, 
for instance, must be limited due to breakdown 
considerations. If the axial electric field becomes too 
large, inter-electrode arcing can occur, and the 
resulting shorts can degrade device performance. 
Furthermore, high heat dissipation in these 
concentrated arc discharges can cause severe erosion of 
wall material and reduction in channel operating life. 

Experience has shown that a realistic limit for the 
axial electric field is about 40 kV/cm. The existence of 
such a design constraint is illustrated by the dashed 
lines marked ±S KUm on the performance map of Fig 5. 
This implies that DCW accelerators operating with 
large Q are confined to low power density operation. 
High power densities are attainable with low Cl at the 
expense of decreased efficiency. 

An additional design constraint arises from the 
observation that excessive heat generation and material 
erosion will occur if the current density entering the 
electrode surface becomes too large. The actual 
limiting value depends on electrode material and 
geometry, but experience indicates good channel 
durability for values up to about 10 A/cm 2 . For high 
power density applications with short-duration 
operating requirements, this value can be exceeded. 

Conclusions 

The classical cross plane averaged performance 
theory for diagonal conducting wall MHD accelerators 
has been extended to include a current dependent 
sheath layer voltage drop. This approach yields 
analytical performance relationships and diagrams that 
can be used to illuminate the rudimentary behavior of 
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these devices and underscore the fundamental interplay 
of basic parameters. This simplified theoretical 
treatment is not intended for detailed performance 
prediction of practical devices but rather to aid in 
developing an intuitive understanding of device 
operation of general value to analysts, designers, and 
experimentalists, alike. Despite its limitations, the 
theory can be extremely useful in defining anticipated 
performance ranges while accounting for critical non- 
ideal effects. 
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